Monitoring tissue temperature while using an irrigated catheter

ABSTRACT

A method, consisting of modeling physical parameters representative of a probe in proximity to body tissue during an ablation procedure performed by the probe. The method also includes measuring a subgroup of the physical parameters during a non-ablation stage of the ablation procedure so as to generate measured non-ablative values of the subgroup, and measuring the subgroup of the physical parameters during an ablation stage of the ablation procedure so as to generate measured ablative values of the subgroup. In response to the modeling, the method includes generating calculated non-ablative values of the subgroup for the non-ablation stage, and generating calculated ablative values of the subgroup for the ablation stage. The method compares the measured non-ablative values with the calculated non-ablative values, and compares the measured ablative values with the calculated ablative values, so as generate optimal values of the physical parameters.

CROSS-REFERENCE TO RELATED APPLICATIONS

This patent application is a continuation of U.S. patent application Ser. No. 15/968,221, filed May 1, 2018, now U.S. Pat. No. 11,103,305, which is a continuation of U.S. patent application Ser. No. 13/036,469, filed Feb. 28, 2011, now U.S. Pat. No. 9,980,772, which claims the benefit of U.S. Provisional Patent Application No. 61/312,447, filed Mar. 10, 2010, each of which is incorporated by reference in its entirety.

FIELD OF THE INVENTION

The present invention relates generally to ablation of body tissue, and specifically to measuring the temperature of the tissue during the ablation.

BACKGROUND OF THE INVENTION

Medical procedures involving ablation of the heart may be used to cure a variety of cardiac arrhythmia, as well as to manage atrial fibrillation. Such procedures are known in the art. Other medical procedures using ablation of body tissue, such as treating varicose veins, are also known in the art. The ablation energy for these procedures may be in the form of radio-frequency (RF) energy.

The application of the ablation energy to body tissue, if uncontrolled, may lead to an unwanted increase of temperature of the tissue. It is consequently important to measure the temperature of the tissue during any medical procedure involving ablation.

SUMMARY OF THE INVENTION

An embodiment of the present invention provides a method, including:

-   -   modeling physical parameters representative of a probe in         proximity to body tissue during an ablation procedure performed         by the probe;     -   measuring a subgroup of the physical parameters during a         non-ablation stage of the ablation procedure so as to generate         measured non-ablative values of the subgroup;     -   measuring the subgroup of the physical parameters during an         ablation stage of the ablation procedure so as to generate         measured ablative values of the subgroup;     -   in response to the modeling:     -   generating calculated non-ablative values of the subgroup for         the non-ablation stage; and     -   generating calculated ablative values of the subgroup for the         ablation stage; and     -   comparing the measured non-ablative values with the calculated         non-ablative values, and comparing the measured ablative values         with the calculated ablative values, so as generate optimal         values of the physical parameters.

Typically, the physical parameters include a matrix of temperatures of the body tissue.

In a disclosed embodiment the subgroup includes a temperature, and the measured non-ablative values include a first time constant of a first change of the temperature, and the measured ablative values include a second time constant of a second change of the temperature.

In another disclosed embodiment the method includes irrigating the tissue via the probe, wherein the subgroup includes an irrigation rate alternating between a first measured irrigation rate during the non-ablation stage and a second measured irrigation rate during the ablation stage.

In yet another disclosed embodiment, the subgroup includes power supplied by the probe, alternating between a first power level able to ablate the tissue during the ablation stage, and a second power level not able to ablate the tissue during the non-ablation stage.

Typically, the subgroup includes at least one of a location, an orientation, and a contact area of the probe with respect to the body tissue.

Typically, the ablation procedure consists of radio-frequency ablation of the body tissue.

There is further provided, according to an alternative embodiment of the invention, a method, including:

-   -   computing a first estimated geometrical vector representative of         a location, an orientation, and a contact area of a probe with         respect to body tissue in response to electrical impedances         between the probe and the body tissue;     -   computing a second estimated geometrical vector of the probe in         response to a set of thermal parameters associated with an         ablation procedure performed on the tissue by the probe;     -   comparing the first and second estimated geometrical vectors to         formulate an optimal geometric vector; and     -   estimating a temperature of the tissue using the optimal         geometric vector.

Typically, the probe includes two or more electrodes insulated from each other, and the electrical impedances are respectively measured between the two or more electrodes and the body tissue.

The first estimated geometrical vector may include a first distance, a first angle, and a first contact area between the probe and the body tissue, the second estimated geometrical vector may include a second distance, a second angle, and a second contact area between the probe and the body tissue, and the optimal geometrical vector may include a third distance, a third angle, and a third contact area between the probe and the body tissue.

In a disclosed embodiment the ablation procedure includes supplying ablation power intermittently to the tissue via the probe, so that the ablation power alternates between a first measured ablation power level sufficient for ablation of the tissue during an ablation stage, and a second measured ablation power level not sufficient for ablation of the tissue during a non-ablation stage.

Typically, computing the second estimated geometrical vector further includes computing the vector in response to the electrical impedances.

Estimating the temperature of the tissue may include estimating a map of temperature values for the tissue.

There is further provided, according to an embodiment of the present invention, apparatus, including:

-   -   a probe, configured to be placed in proximity to body tissue;         and     -   a controller which is configured to:     -   compute a first estimated geometrical vector representative of         an orientation, and a contact area of the probe with respect to         the body tissue in response to electrical impedances between the         probe and the body tissue,     -   compute a second estimated geometrical vector of the probe in         response to a set of thermal parameters associated with an         ablation procedure performed on the tissue by the probe,     -   compare the first and second estimated geometrical vectors to         formulate an optimal geometric vector, and     -   estimate a temperature of the tissue using the optimal geometric         vector.

There is further provided, according to an alternative embodiment of the present invention, apparatus, including:

-   -   a probe, configured to be placed in proximity to body tissue;         and     -   a controller which is configured to:     -   model physical parameters representative of the probe during an         ablation procedure performed by the probe on the body tissue;     -   measure a subgroup of the physical parameters during a         non-ablation stage of the ablation procedure so as to generate         measured non-ablative values of the subgroup;     -   measure the subgroup of the physical parameters during an         ablation stage of the ablation procedure so as to generate         measured ablative values of the subgroup;     -   in response to the modeling:     -   generate calculated non-ablative values of the subgroup for the         non-ablation stage; and     -   generate calculated ablative values of the subgroup for the         ablation stage; and     -   compare the measured non-ablative values with the calculated         non-ablative values, and compare the measured ablative values         with the calculated ablative values, so as generate optimal         values of the physical parameters.

There is further provided, according to an embodiment of the present invention, a computer software product consisting of a non-transitory computer-readable medium having computer program instructions recorded therein, which instructions, when read by a computer, cause the computer to:

-   -   model physical parameters representative of a probe in proximity         to body tissue during an ablation procedure performed by the         probe;     -   measure a subgroup of the physical parameters during a         non-ablation stage of the ablation procedure so as to generate         measured non-ablative values of the subgroup;     -   measure the subgroup of the physical parameters during an         ablation stage of the ablation procedure so as to generate         measured ablative values of the subgroup;     -   in response to the modeling:     -   generate calculated non-ablative values of the subgroup for the         non-ablation stage; and     -   generate calculated ablative values of the subgroup for the         ablation stage; and     -   compare the measured non-ablative values with the calculated         non-ablative values, and compare the measured ablative values         with the calculated ablative values, so as to generate optimal         values of the physical parameters.

The present disclosure will be more fully understood from the following detailed description of the embodiments thereof, taken together with the drawings, in which:

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic, pictorial illustration of a catheter probe ablating system, according to an embodiment of the present invention;

FIG. 2 is a schematic cross-section of a distal tip of a catheter probe used in the system, according to an embodiment of the present invention;

FIGS. 3 and 4 are schematic theoretical diagrams illustrating parameters used by the ablating system, according to embodiments of the present invention;

FIG. 5 is a schematic theoretical diagram illustrating a method used by the ablating system to estimate values for a geometric vector, according to an embodiment of the present invention;

FIG. 6 is a schematic block diagram illustrating a method of analysis, according to an embodiment of the present invention;

FIG. 7 illustrates schematic graphs generated by the ablating system, according to an embodiment of the present invention; and

FIG. 8 is a flowchart showing steps followed by a controller of the ablating system, according to an embodiment of the present invention.

DETAILED DESCRIPTION OF EMBODIMENTS Overview

An embodiment of the present invention uses a mathematical model of physical parameters related to an ablation procedure. The modeled parameters represent the distal tip of a probe in proximity to body tissue during the procedure, which is performed by the probe. The ablation procedure consists of an ablation stage, in which power is applied to the body tissue at a level sufficient to ablate the tissue, and a non-ablation stage in which the power level is reduced so that no tissue ablation occurs. Typically, the period of the non-ablation stage is brief in comparison with the period of the ablation stage. During the ablation stage, a subgroup of the physical parameters is measured, generating ablative values of the subgroup. During the non-ablation stage, the subgroup of parameters is measured to generate non-ablative values of the subgroup.

The model is used to calculate ablative values of the subgroup for the ablation stage, and non-ablative values for the non-ablation stage. The measured and calculated values for the two stages are compared, so as to generate optimal values of the physical parameters.

The physical parameters that are modeled include values of power supplied by the probe, for ablating the tissue. The power may be in substantially any form, such as radio-frequency (RF), optical, or ultrasonic power. The parameters also include, but are not limited to, parameters such as a contact force of the probe on the body tissue, an irrigation rate of cooling fluid provided by the probe, impedances/conductivities of the tissue, the cooling fluid and blood surrounding the tissue, temperatures of elements in proximity to the probe, including temperatures at the distal tip, and geometrical parameters of the probe.

The model may be used to evaluate ablation of heart tissue. In one embodiment a probe having multiple electrodes is inserted into proximity with the tissue, and one or more of the electrodes transfers RF power to the tissue in the ablation procedure described above. During the procedure impedances between the electrodes and the tissue are measured, and the impedances are used to estimate a location, an orientation, and/or a contact area of the probe's distal tip with respect to the tissue. Typically, irrigation rates of the cooling fluid are altered in synchrony with the ablation procedure periods of ablation and non-ablation.

The estimate may be provided to the model which uses the estimate, as well as values of other parameters (referred to above) that are measured directly or indirectly during the ablation procedure, to produce a temperature map of the tissue being ablated.

In addition to producing the temperature map, the model may be configured to revise the estimate of the distal tip's location, orientation, and/or contact area. Typically, the process of generating and revising the estimate is performed on an iterative basis until optimal values of the location, orientation and/or contact area are attained. The model may then use the optimal values to provide a final formulation of the temperature map of the tissue.

A similar iterative process to that described above (for the probe's location, orientation, and/or contact area) may be applied for other parameters of the model.

System Description

Reference is now made to FIG. 1 , which is a schematic, pictorial illustration of a catheter probe ablating system 20, and to FIG. 2 which is a schematic cross-section of a distal tip of a catheter probe 22 used in the system, according to embodiments of the present invention. In system 20, probe 22 is inserted into a lumen 23, such as a chamber of a heart 24, of a subject 26. The probe is used by an operator 28 of system 20, during a procedure which includes performing ablation of body tissue 25.

The functioning of system 20 is managed by a system controller 30, comprising a processing unit 32 communicating with a memory 34, wherein is stored software for operation of system 20. Controller 30 is typically an industry-standard personal computer comprising a general-purpose computer processing unit. However, in some embodiments, at least some of the functions of the controller are performed using custom-designed hardware and software, such as an application specific integrated circuit (ASIC) or a field programmable gate array (FPGA). Controller 30 is typically operated by operator 28 using a pointing device 36 and a graphic user interface (GUI) 38, which enable the operator to set parameters of system 20. GUI 38 typically also displays results of the procedure to the operator.

The software in memory 34 may be downloaded to the controller in electronic form, over a network, for example. Alternatively or additionally, the software may be provided on non-transitory tangible media, such as optical, magnetic, or electronic storage media.

FIG. 2 illustrates in cross-section a distal tip 40 of probe 22 in proximity to tissue 25. Distal tip 40 is assumed to be immersed in fluid 27, so that tissue 25 has a surface 29 contacting the fluid. (As described below, fluid 27 typically comprises a mixture of blood and saline solution.) By way of example, distal tip 40 is assumed herein to be formed from an insulating substrate 41 in the shape of a cylinder 45 closed by a generally flat surface 47 at one end. Cylinder 45 has an axis of symmetry 51. As shown in FIG. 2 , a curved section 49 joins flat surface 47 and cylinder 45. A typical diameter of cylinder 45 is 2.5 mm, and a typical radius of the curved section is 0.5 mm.

Distal tip 40 comprises a first electrode 44, a second electrode 46, and a third electrode 48, the electrodes being insulated from each other. The electrodes typically comprise thin metal layers formed over insulating substrate 41. Typically, the distal tip has other electrodes, insulated from electrodes 44, 46, and 48, which for simplicity are not shown in the diagram. Tip electrode 44 has the shape of a cup with a flat base, and is herein also referred to as the cup electrode. Cup electrode 44 typically has a thickness in a range from approximately 0.1 mm to approximately 0.2 mm. Second and third electrodes 46 and 48 are usually in the form of rings, and are also referred to herein as ring electrodes 46 and 48. In the present disclosure, electrodes 44, 46, and 48, and other electrodes of the distal tip, are also referred to herein collectively as electrodes 53.

Electrodes 44, 46, and 48 are connected to system controller 30 by wires, not shown in the figures. As described below, at least one of the electrodes is used to ablate tissue 25. An ablation module 50, under control of controller 30, is configured to be able to set and measure a level of the ablation power delivered by each of the electrodes.

Typically, during the ablation, heat is generated in the ablating electrode and in the surrounding region. In order to dissipate the heat, system 20 irrigates the region, using a number of small irrigation apertures 52 in the cup electrode. The apertures typically have diameters in an approximate range 0.1-0.2 mm. An irrigation tube 54 supplies saline solution to the apertures, and the rate of flow of the saline solution through the apertures (causing fluid 27 to be a mixture of blood and saline solution) is controlled by an irrigation module 56 in the system controller. The saline rate of flow is typically in the range of approximately 10-20 cc/minute, but may be higher or lower than this range.

A saline temperature sensor 58, typically a thermocouple, is located in tube 54, and provides a signal to module 56 allowing the module to measure a temperature T_(S) of the saline solution input to apertures 52. While the saline solution may be provided by module 56 at room ambient temperature, e.g., in a range of approximately 19-22° C., the solution is heated during its flow in probe 22, so that temperature T_(S) is typically in an approximate range of 26-28° C.

Typically, one or more location sensing devices 61 are incorporated in the distal tip. Devices 61 are configured to provide signals to controller 30 enabling the unit to ascertain the position and/or orientation of distal tip 40, The position and orientation are typically measured with respect to subject 26. One such location sensing device comprises the magnetic navigation CARTO system, produced by Biosense-Webster Inc, of Diamond Bar, CA.

Alternatively or additionally, the location sensing device operates by measuring currents between electrodes 53 of the distal tip and one or more electrodes on the skin of subject 26. A skin electrode 62 is assumed to be attached to subject 26. Such a method of measuring the location of distal tip 40 is disclosed in U.S. Patent Application 2010/0079158, which is incorporated herein by reference. Thus, in addition to serving as electrodes providing ablation, electrodes 53 may also be used for other functions, such as for electrophysiological sensing and/or for locating the distal tip.

As necessary, when used for other functions, controller 30 may differentiate between the currents for the different functions by frequency multiplexing. For example, ablation module 50 typically generates its ablation power at frequencies of the order of hundreds of kHz, while the position sensing frequencies may be at frequencies of the order of 1 kHz.

Typically, distal tip 40 comprises other elements, which for clarity and simplicity are not shown in FIG. 2 . One such element comprises a force sensor, which is configured to measure a force exerted by tip 40 on tissue 42.

Distal tip 40 comprises one or more generally similar temperature sensors 82 which are fixedly connected, by an insulator, to the outer surface of cup electrode 44, so as to protrude from the surface. Sensors 82 have a typical diameter of approximately 0.3 mm and a length of approximately 1.5 mm. In one embodiment sensors 82 are thermistors NTC Type AB6, produced by General Electric Company of Schenectady, New York. By way of example, the following description assumes there are three sensors 82 (only two are shown in the diagram) symmetrically distributed with respect to axis 51, and located on a curved section 86 of the cup electrode. Curved section 86 of the cup electrode overlays curved section 49 of the distal tip.

A magnified section 88 of FIG. 2 illustrates one of sensors 82 in more detail. As shown in section 88, an insulator 84 separates sensor 82 from curved section 86 of the cup electrode. Insulator 84 is selected to provide good thermal and electrical insulation, and in some embodiments insulator 84 may comprise an adhesive that bonds sensor 82 to curved section 86. Wires 90 connect sensor 82 to controller 30.

By having sensors 82 protrude from the outer surface of cup electrode 44, the sensors are able to intimately contact tissue 25. Controller 30 is thus able to use signals from the sensors to provide direct temperature measurements of the tissue.

FIGS. 3 and 4 are schematic theoretical diagrams illustrating parameters used by system 20, according to an embodiment of the present invention. FIG. 3 illustrates a situation when tip 40 is above surface 29 of tissue 25, and FIG. 4 illustrates a situation when tip 40 is in contact with the tissue. System 20 estimates the temperature of tissue 25 during an ablation procedure, by controller 30 measuring and analyzing values of thermal parameters, illustrated in the figures, that are associated with the procedure.

The analysis considers the rate of energy transfer in a small volume 100 of fluid 27 surrounding cup electrode 44 and bordered by surface 29 of tissue 25. The analysis also considers the rate of energy transfer in a small volume 102 in the tissue adjacent to volume 100.

In general, an equation relating areas of distal tip 40 is: A(t)=A(b)+A(c)  (1)

-   -   where A(t) is an overall area of distal tip 40 in contact with         volume 100,     -   A(b) is the part of area A(t) that is exposed to the fluid, and     -   A(c) is a contact area between the tip and the tissue.

In the situation illustrated in FIG. 3 , A(c)=0, so that A(b)=A(t). In the situation illustrated in FIG. 4 , A(c)>0, so that A(b)<A(t).

The power into volume 100 is assumed to be caused by electromagnetic radio-frequency (RF) ablation power Q_(RF-L) dissipated in the liquid (the blood and the saline solution) and by heat power Q_(T) transferred into volume 100 from volume 102 of tissue 25, due to the higher temperature of the tissue compared to the temperature of blood, T_(B). (Typically, the higher temperature of the tissue may reach approximately 90-100° C.)

The value Q_(RF-L) of the ablation power dissipated in the blood and saline solution of volume 100 is a function of the distal tip area exposed to the fluid, A(b), the conductivity of the blood, σ_(B)(T), and of the conductivity of the saline solution, σ_(S)(T). (Both conductivities are strongly dependent on the temperature T since both have electrolytic properties. Typically the conductivities change by about 2%/° C.)

Q_(RF-L) may be written: Q _(RF_L) =Q _(RF)(σ_(B)(T),σ_(S)(T)),A(b))  (2)

The ablation power dissipated in volume 102 (discussed in more detail below) of the tissue causes the tissue to rise in temperature above a standard body temperature assumed to be 37° C. Thus, there is an array, or map, of tissue temperatures in volume 102 of the tissue. The map of tissue temperatures may be represented by a matrix [T], where each element of the matrix is a temperature of a respective volume element in the region. The heat power transferred from volume 102 to volume 100, Q_(T), is a function of the map of temperatures [T] and T_(B).

Matrix [T] is illustrated schematically in FIG. 4 , wherein a line 104 represents a 100° C. isotherm, and a line 106 represents a 50° C. isotherm. Regions with temperatures of 50° C. or more typically correspond to necrotic lesions.

A power, P_(in), into volume 100 is given by an equation (3): P _(in) =Q _(RF)(σ_(B)(T),σ_(S)(T),A(b))+Q _(T)([T _(T) ],T _(B))  (3)

where Q_(T)([T_(T)],T_(B)) is a function representing the heat power transferred into volume 100 due to differences between blood temperature T_(B) and the temperatures of [T_(T)], typically comprising the temperatures of the hottest part of the tissue such as those within line 104.

The power out of volume 100 is carried by the blood mass flow and by the irrigating saline solution mass flow. An expression for the power dissipated by the fluids, P_(out), is:

$\begin{matrix} {P_{out} = {{\frac{{dm}_{B}}{ac} \cdot C_{pB} \cdot \left( {T_{out} - T_{B}} \right)} + {X \cdot \frac{{dm}_{S}}{ac} \cdot C_{ps} \cdot \left( {T_{out} - T_{s}} \right)}}} & (4) \end{matrix}$

where

$\frac{{dm}_{B}}{dt}$ is the blood mass flow,

-   -   C_(pB) is the specific heat capacity of the blood,

$\frac{{dm}_{s}}{dt}$

-   -    is the saline mass flow,     -   C_(pS) is the specific heat capacity of the saline solution,     -   T_(S) is the input temperature of the saline solution,     -   T_(B) is the temperature of the blood entering volume 100,         typically assumed to be constant at 37° C.,     -   T_(out) is the temperature of the blood and saline solution         mixture leaving volume 100, and     -   X is a factor depending on the amount of mixing of the saline         solution and the blood in volume 100.

Typically, temperature T_(out) may be determined by sensors 82. In the situation illustrated by FIG. 3 , T_(out) may usually be measured using all three sensors; in the situation illustrated by FIG. 4 , and assuming an angle θ (defined below) is greater than approximately 45°, T_(out) is typically registered by at least one of sensors 82.

Factor X in equation (4) is a function of whether the saline solution exits from irrigation apertures 52 as a turbulent or a laminar flow, since the type of flow governs the length of the saline jet exiting from each aperture. The type of flow in turn is dependent on a diameter d of the apertures and the saline mass rate of flow

$\frac{{dm}_{s}}{dt}.$ The factor X is also dependent on the proximity of distal tip 40 to surface 29 of tissue 25 i.e., on a distance Δx between the distal tip and the surface. As illustrated, the distance Δx may be positive or negative, depending whether the distal tip is above surface 29 (FIG. 3 ), or below the surface (FIG. 4 ).

The factor X is further dependent on the directions of the flow from the apertures relative to the tissue. The directions of flow are a function of an angle θ made by distal tip 40 with tissue 25.

Thus, an energy balance equation (4) may be written:

$\begin{matrix} {P_{in} = {P_{out} = {{\frac{{dm}_{B}}{dt} \cdot C_{pB} \cdot \left( {T_{out} - T_{B}} \right)} + {X{\frac{{dm}_{s}}{dt} \cdot C_{ps} \cdot \left( {T_{out} - T_{s}} \right)}}}}} & (5) \end{matrix}$

Considering volume 102, an ablation power Q_(RF-TIS) dissipated in the tissue may be expressed as the equation: Q _(RF_TIS)=Σ_(V) q _(TIS)(σ_(TIS))  (6)

where each volume element V of volume 102 has an RF power q_(TIS) dissipated in the element, according to a potential across the element and a current through the element, and where each volume element has a respective conductance σ_(TIS).

The two ablation powers, Q_(RF-L) and Q_(RF-TIS), are supplied by ablation module 50, and are related by the equation: Q _(RF) =Q _(RF-L) +Q _(RF-TIS)  (7)

Q_(RF) is both controlled and measured by controller 30. Typically, during an ablation procedure, Q_(RF) is of the order of 30 W.

As is described further below, embodiments of the present invention evaluate the expressions of equations (1)-(7). The evaluations may be used to generate values for elements of matrix [T]. In addition, the evaluations may be used to estimate values for Δx, θ, and/or A(c). In the description herein, the ordered triple (Δx, θ, A(c)) is referred to as the geometrical vector of the distal tip.

FIG. 5 is a schematic theoretical diagram illustrating an alternative method used by system 20 to estimate values for Δx, θ, and/or A(c), according to an embodiment of the present invention. For simplicity, only elements of system 20 that are required for explaining the alternative method are labeled in the figure. In the alternative method, controller 30 measures an impedance X_(C) between cup electrode 44 and skin electrode 62 (of subject 26). In addition, the controller measures an impedance X_(R) between ring electrode 48 and the skin electrode.

The impedances X_(C), X_(R), are functions of the conductivity of fluid 27 in which both electrodes are immersed, and of their different effective path lengths in the fluid to surface 29. The conductivity of fluid 27 is dependent on the blood and saline conductivities, σ_(B)(T), σ_(S)(T). Both impedances are also functions of the conductivity of the subject tissue between surface 29 and skin electrode 62, which have substantially the same path lengths for both impedances.

As stated above, distal tip 40 typically comprises more than two electrodes, so that the explanation above may be applied for impedances between each electrode 53 and skin electrode 62.

As described below, controller 30 uses the measured impedances (after normalizing to take into account the different surface areas of the cup electrode and the one or more ring electrodes) to estimate values for Δx, θ, and/or A(c) for distal tip 40. The controller typically uses, inter alia, values of other physical parameters, such as conductivities of tissue, blood, and saline solution, to make the estimation.

FIG. 6 is a schematic block diagram 200 illustrating a method of analysis used by controller 30, according to an embodiment of the present invention.

A controllable inputs block 202 represents the inputs to distal tip 40 which controller 30 is configured to direct, and which the unit is able to measure. As explained above, these inputs comprise the overall ablation power generated in ablation module 50, and dissipated in volumes 100 and 102 (FIGS. 3 and 4 ), and the saline irrigation flow into distal tip 40.

A dimensions block 204 represents the mechanical dimensions of distal tip 40. Thus, assuming probe 22 is cylindrical, the dimensions comprise a diameter of the cylinder and a length of the distal tip that is used by the controller in its analysis.

The analysis by the controller uses values of physical parameters that are known in the art, and that may be saved in memory 34 for use by controller 30. A physical parameters block 206 represents these parameters, which comprise, for example, tissue and blood electrical conductivities, and rates of change of these conductivities with temperature. The parameters also comprise tissue thermal conductivity and its rate of change with temperature.

During an ablation procedure controller 30 is able to measure temperatures provided by saline sensor 58 and distal tip sensors 82. Typically at least one of sensors 82 is in direct contact with tissue 25 and directly measures a tissue temperature a few tenths of a millimeter below surface 29. The controller is also able to measure the impedance presented to ablation module 50 during the procedure. In addition, the time of each of these measurements may be registered by the controller, and used to determine rates of temperature and impedance changes. A sensors block 208 represents thermal and non-thermal parameter measurements, such as those described above, which may be registered by the controller.

In a modeling process, represented by a model processing block 210, the controller incorporates the inputs of blocks 202, 204, 206, and 208 into equations (1)-(7). The modeling process typically applies finite element methods (FEMs) and computational fluid dynamics (CFD), in order to estimate values of unknown parameters of the equations, as well as to generate results of the equations. Both FEM and CFD are known in the art, and the application of these processes to analyze equations (1)-(7) will be apparent to one having ordinary skill in the art.

The controller evaluates the results of the equations, and from an analysis of the results the unit may make an estimate of thermally-dependent parameters, herein also termed thermal parameters, affecting a geometric estimation block 214, described below. The thermal parameters, represented by a model parameters block 212, typically include evaluations of the blood and saline conductivities σ_(B)(T), σ_(S)(T), from estimates of temperature T. From the analysis the controller may also make an estimate of the map of temperatures [T], as represented by a temperature map block 215.

In addition to the modeling process represented by block 210, controller 30 implements a geometric estimation of vector (Δx, θ, A(c)), using the process described above with reference to FIG. 5 . The estimation is represented by geometry estimation block 214. Block 214 receives as inputs, represented by an impedance value block 213, the impedances (measured by controller 30) between electrodes 53 and skin electrode 62. While the estimation represented by block 214 may be independently made by the controller, i.e., using only the values of block 213, the controller may apply some of the thermal parameters of block 212 in performing the estimation, as indicated by an arrow 216.

Controller 30 transfers the output of block 214, represented by a geometry output block 218, to the modeling process of block 210. Typically, as described below with reference to flowchart 300, the controller iterates between blocks 210 and 218 to generate an optimal geometric vector (Δx_(O), θ_(O), A(c)_(O)). Thus, while an initial estimate of the vector may be substantially independent of the thermal parameters of block 212, a subsequent estimate uses the thermal parameters.

FIG. 7 illustrates schematic graphs generated by system 20, according to an embodiment of the present invention. The graphs depict simulated measurements during an intermittent ablation procedure for ablating tissue. The ablation procedure comprises ablation stages interleaved with non-ablation stages. Graphs 250 and 252 are temperature vs. time graphs, a graph 254 is of the rate of saline solution flow vs. time, and a graph 256 is of the ablation power vs. time.

During the procedure the controller applies a known ablation power, 30 W, and a known irrigation flow rate 18 cc/min, between an initial time (0 s) and a time T1, for a period ΔT₁. Period ΔT₁ corresponds to an ablation stage of the procedure, since the power supplied during period ΔT₁ is sufficient to cause necrosis of the tissue, i.e., is sufficient to ablate the tissue. At time T1, the controller switches off the ablation power and reduces the irrigation flow to 2 cc/min for a short time period ΔT₂, which finishes at a time T2. In one embodiment ΔT₂ is typically approximately 500 ms. During period ΔT₂, corresponding to a non-ablation stage of the procedure, the tissue does not ablate. At time T2 the controller applies the ablation powers and irrigation rates used during period ΔT₁, until time T3, for a period ΔT₃, which is an ablation stage. From time T3, the ablation power is switched off, and the irrigation rate is reduced to the same rate as for period ΔT₂. The period from T3 is a non-ablation stage.

Graph 252 shows the temperature measured by a sensor 82 of tip 40, which by way of example is assumed to be in contact with surface 29 of tissue 25. During time periods ΔT₁ and ΔT₃, i.e., during ablation stages wherein the tissue is being ablated, the controller is able to extract from the temperature and time measurements corresponding time constants τ₁, τ₃, indicative of the positive rates of change of temperature during these periods. During time period ΔT₂, i.e., during non-ablation stages wherein the tissue is not being ablated, the controller evaluates a time constant τ₂, indicative of the negative rate of change of temperature during period ΔT₂. The controller may also evaluate a time constant τ₄ for the non-ablation stage following time T3.

The controller may use the values of the time constants in the analysis performed in modeling block 210 in order to determine a corresponding temperature map [T] of the tissue being ablated. Graph 250 shows the temperature of a hot region of the tissue determined by modeling block 210.

During ablation of subject 26, the controller may perform the intermittent, “on/off,” ablation procedure illustrated by graphs 254 and 256. During such an ablation procedure, the controller may measure the temperature near distal tip 40 with sensors 82, and use the temperature measurements to derive time constants for the temperature change while ablation is occurring and while it is not being performed, as exemplified above. The time constants may be used in modeling block 210 to derive a temperature map for the tissue of subject 26.

The intermittent on/off ablation procedure described above has assumed a specific value of the ablation power and a specific value of irrigation flow during the “on” stage of the procedure. Similarly, the procedure assumed another specific value of the ablation power, in this case zero, and another specific value of irrigation flow during the “off” stage of the procedure. It will be understood that these values are by way of example, and that other power and flow rate values may be implemented, and measured, by controller 30 during an on/off ablation procedure. Similarly, controller 30 may configure other values than those exemplified above for the ablation stage periods and for the non-ablation stage periods of the procedure. Using different values may allow the controller, in the analysis performed in block 210, to produce a temperature map [T} with increased accuracy.

FIG. 8 is a flowchart 300 showing steps followed by controller 30, according to an embodiment of the present invention. In an initial parameter storing step 302, the unit stores values corresponding to dimensions block 204 and physical parameters block 206 (FIG. 6 ) in memory 34.

In an initiate procedure step 304, operator 28 inserts probe 22 into heart 24, and positions distal tip 40 in a desired location of the heart using location devices 61.

In an operational step 306, the operator begins an ablation procedure. During the procedure controller 30 adjusts and measures the ablation power provided by module 50, as well as the rate of saline irrigation provided by module 56. The procedure typically involves the intermittent on/off procedure described above with reference to FIG. 7 . The parameters of controllable inputs of block 202, comprising parameters of the on and off states of the on/off procedure, may be set by operator 28. The parameters of the on and off states comprise the time periods of each state, as well as the delivered power level and saline flow rate during the states.

Alternatively or additionally, controller 30 may be configured to set the parameters of the controllable inputs of block 202. The controller may set the parameters semi-automatically, i.e. with some input from the operator, or automatically, i.e., with substantially no input from the operator. In some embodiments, the controller may set the parameters adaptively, typically in response to results derived from modeling block 210.

Also in step 306, controller 30 records the sensor input values, corresponding to sensor measurements block 208.

In an impedance measuring step 308, which is typically implemented concurrently with step 306, the controller measures impedances, described above with reference to FIG. 5 , between electrodes 53 and skin electrode 62. In the case of the intermittent ablation procedure described above, the impedances may be measured during the on state and during the off state.

In a vector estimation step 310, controller 30 uses the impedances of step 308 to estimate a value for the probe geometric vector (Δx, θ, A(c)). The estimation typically uses physical parameters, such as conductivities, corresponding to physical parameters block 206. In some embodiments, the controller may use the impedances to measure values of a subgroup of the elements of the geometric vector.

In a model application step 312 the values estimated in step 310 are transferred to model processing block 210, and the controller applies the processing of block 210 to make a further estimate of the geometric vector.

In a comparison step 314 the controller compares the estimates of steps 310 and 312. If the difference between the two estimates is greater than a pre-defined limit, the flowchart returns to step 306. If the difference is less than or equal to the pre-defined limit, the flowchart continues to a final step 316.

In final step 316, the controller formulates an optimal geometric vector (Δx_(O), θ_(O), A(c)_(O)), from the vectors determined in steps 310 and 312. The controller uses the optimal geometric vector, in model processing block 210, to generate temperature map [T], corresponding to temperature map block 215. The controller typically uses a similar iterative process to that described for generating the optimal geometric vector, by comparing physical measurements of sensors such as sensors 82 with modeled values, in order to optimize temperature map [T].

While the description above assumes that the ablation power is typically RF power, it will be understood that embodiments of the present invention may use substantially any type of power for ablation, such as, but not limited to, optical or ultrasonic power.

It will thus be appreciated that the embodiments described above are cited by way of example, and that the present invention is not limited to what has been particularly shown and described hereinabove. Rather, the scope of the present invention includes both combinations and subcombinations of the various features described hereinabove, as well as variations and modifications thereof which would occur to persons skilled in the art upon reading the foregoing description and which are not disclosed in the prior art. 

We claim:
 1. A computer software product comprising a non-transitory computer-readable medium having computer program instructions recorded therein, which instructions, when read by a computer, cause the computer to: model physical parameters representative of a probe in proximity to body tissue during an ablation procedure performed by the probe on the body tissue by generating calculated non-ablative values of a subgroup for the non-ablation stage and generating calculated ablative values of the subgroup for the ablation stage; measure a subgroup of the physical parameters through at least one of the sensor and the electrode during the non-ablation stage of the ablation procedure so as to generate measured non-ablative values of the subgroup; measure the subgroup of the physical parameters through at least one of the sensor and the electrode during the ablation stage of the ablation procedure so as to generate measured ablative values of the subgroup; and compare the measured non-ablative values with the calculated non-ablative values, and compare the measured ablative values with the calculated ablative values, so as to generate optimal values of the physical parameters.
 2. The computer software product according to claim 1, wherein the physical parameters comprise a matrix of temperatures of the body tissue.
 3. The computer software product according to claim 1, wherein the subgroup comprises a temperature, and wherein the measured non-ablative values comprise a first time constant of a first change of the temperature, and wherein the measured ablative values comprise a second time constant of a second change of the temperature.
 4. The computer software product according to claim 1, and comprising irrigating the tissue via the probe, and wherein the subgroup comprises an irrigation rate alternating between a first measured irrigation rate during the non-ablation stage and a second measured irrigation rate during the ablation stage.
 5. The computer software product according to claim 1, wherein the subgroup comprises power supplied by the probe, alternating between a first power level able to ablate the tissue during the ablation stage, and a second power level not able to ablate the tissue during the non-ablation stage.
 6. The computer software product according to claim 1, wherein the subgroup comprises at least one of a location, an orientation, and a contact area of the probe with respect to the body tissue.
 7. The computer software product according to claim 1, wherein the ablation procedure comprises radio-frequency ablation of the body tissue.
 8. A method, comprising: modeling physical parameters representative of a probe in proximity to body tissue during an ablation procedure performed by the probe; measuring a subgroup of the physical parameters during a non-ablation stage of the ablation procedure so as to generate measured non-ablative values of the subgroup; measuring the subgroup of the physical parameters during an ablation stage of the ablation procedure so as to generate measured ablative values of the subgroup; in response to the modeling: generating calculated non-ablative values of the subgroup for the non-ablation stage; and generating calculated ablative values of the subgroup for the ablation stage; and comparing the measured non-ablative values with the calculated non-ablative values, and comparing the measured ablative values with the calculated ablative values, so as generate optimal values of the physical parameters.
 9. The method according to claim 8, wherein the physical parameters comprise a matrix of temperatures of the body tissue.
 10. The method according to claim 8, wherein the subgroup comprises a temperature, and wherein the measured non-ablative values comprise a first time constant of a first change of the temperature, and wherein the measured ablative values comprise a second time constant of a second change of the temperature.
 11. The method according to claim 8, and comprising irrigating the tissue via the probe, and wherein the subgroup comprises an irrigation rate alternating between a first measured irrigation rate during the non-ablation stage and a second measured irrigation rate during the ablation stage.
 12. The method according to claim 8, wherein the subgroup comprises power supplied by the probe, alternating between a first power level able to ablate the tissue during the ablation stage, and a second power level not able to ablate the tissue during the non-ablation stage.
 13. The method according to claim 8, wherein the subgroup comprises at least one of a location, an orientation, and a contact area of the probe with respect to the body tissue.
 14. The method according to claim 8, wherein the ablation procedure comprises radio-frequency ablation of the body tissue. 